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The  following  report  is  a  summary  of  work  performed  under  the 
above  contract  relating  to  me  ana  of  obtaining  higher  power  output  from 
CO^  lasers. 

A»  Doppler  Broadening 

In  an  effort  to  study  the  effect  of  higher  order  transverse  and 
longitudinal  mode  excitation  on  single  mode  power  output  from  a  CO^ 
laser.  Lax,  Louis  ell  and  McKnight^  [LLM]  have  shown  that  for  a 
homogeneously  broadened  laser  the  fundamental  mode,  through  non - 
linearities,  will  excite  higher  order  modes.  This  means  that  atoms 
which  would  normally  supply  energy  to  the  fundamental  mode  are 
wasted.  It  was  therefore  thought  to  be  of  interest  to  study  the  effect 
of  Doppler  broadening  to  see  if  the  molecular  motion  could  aid  in 
suppressing  higher  order  modes. 

The  problem  of  solving  the  coupled  equations  of  motion  for  the 
atomic  inversion,  atomic  polarization,  and  the  amplitudes  of  the 
various  excited  transverse  modes  for  atoms  having  a  Maxwell -Boltzman 
velocity  distribution  is  very  difficult.  The  usual  technique  of  Fourier 
analyzing  the  variables  in  time  harmonics  works  if  the  atoms  are 
either  1)  motionless,  or  2)  moving  through  a  spatially  harmonic  field 
{i.e. ,  moving  in  the  z  direction  through  purely  longitudinal  modes). 

The  technique  does  not  work  for  atoms  moving  through  the  Hermits - 


Gaussian  fields  of  the 


Jfe¥#:tjra^veree  modes  of  spherical  resonators, 

r:e: 

w-'  -  /’.r _ 


(\)  .M ,  Lax,  W.  K.  Louioell  and  Yf.B.  McKnight,  Jour.  App.  Phys., 
41  3136  (19721.  , 

6 


The  lasers  of  interest  at  Redstone  Arsenal,  however,  are  of 
fairly  high  pressure  and  it  is  easy  to  show  that  for  pressures  on  the 
order  of  70  torr  {1C*  dynes  fern'),  the  mean  free  path  for  any  molecule 
is  small  compared  to  the  wavelength  of  the  laser  radiation  (10.  6}i ). 

The  mean  free  path  is 

la.  1  -  kT  . 

where  o  =  rrd^  is  the  cross  section  for  scattering.  At  room  tempera  - 
o 

ture  and  76  lorr,  we  have 

4  „  M  r  ’o:16U,3  x.,102)  „  2 -4  x  10-4  cm. 

Jz  x  105  x  rr  x  (2  x  10”  )** 

The  wavelength  of  the  laser  radiation  is 

X  =  10.6|js?10  x  10”*cm. 

Therefore,  at  a  pressure  of  76  torr,  the  CO^  molecules  cannot  move 
more  than  a  wavelength  before  a  phase -interrupting  collision  occurs. 
Hence,  we  may  use  the  theory  of  Lax,  Louisell  and  McKright  for 
homogeneously  broadened  lasers  for  those  CO^  lasers  of  interest 
at  Redstone  Arsenal  since  Doppler  broadening  will  not  assist  in  getting 
higher  power  output  by  supplying  additional  atoms  to  the  fundamental 
mode  region. 

The  homogeneous  theory  of  LLM  gave  a  criterion  for  suppressing 
higher  order  transverse  modes.  It  showed  that  large  fundamental 
mode  volumes  are  desirable.  Siegman  and  coworkers  have  been 
studying  diverging  mirrors  to  achieve  this  goal.  However,  these  lead 
to  undesirably  large  diffraction  losses.  As  a  compromise,  we  have 
studied  mirrors  which  have  diverging  spherical  mirrors  in  the  central 
region  of  the  beam  surrounded  by  converging  spherical  mirrors.  We 
expected  that  the  fundamental  mode  volume  would  be  increased  by  the 
diverging  central  part  of  the  mirrois  while  ihe  diffraction  losses 

r '  a 


would  be  reduced  by  the  converging  rime  of  the  mirrors.  The 
following  section  gives  t..s,  -?sult8  of  studies. 

B.  Rcsona'.or  Problem 

The  work  on  the  resor.at.or  problem  is  now  complete  and  the 
results  may  be  found  in  the  enclosed  paper.  The  paper  will  be 
sumbitted  shortly  for  publication.  A  talk  was  presented  at  the 
International  Quantum  Electronics  Meeting  in  Montreal  in  May.  The 
computer  program  now  works  well  and  a  copy  of  the  program  has 
been  sent  to  Dr.  VV.  McKnight  at  Redstone  Arsenal  for  any  further 
geometries  that  need  be  studied.  Dr.  McKnight  is  in  the  process  of 
checking  the  theory  experimentally. 

C.  Two  Mode  Problem 

To  obtain  a  feel  for  how  the  power  output  of  a  multimode  laser 
might  compare  to  that  of  a  single  mode  laser,  we  have  attempted  a 
study  of  a  two  mode  laser.  (Both  modes  are  transverse  modes.  )  Wc 
have  been  rathex  busy  preparing  the  resonator  paper  for  publication 
so  that  not  too  much  progress  has  been  made  on  the  two  mode  problem. 
Wo  have,  however,  been  able  to  derive  the  coupled  equations  of  motion 
for  the  amplitudes  of  the  two  modes,  but  we  have  essentially  been 
limited  by  the  complexity  of  the  integrals  obtained. 

i*or  the  analysis  of  the  two  mode  problem,  we  have  proceeded 
along  the  lines  of  the  (LLM)  paper,  only  we  have  included  two  modes 
rather  tha°  one.  Wo  begin  with  their  coupled  equations  of  motion  for 
the  amplitude  of  any  mode  i,  and  the  inversion  in  the  laser  medium, 
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Here,  b^ ,  us^,  and  are  the  amplitude,  cold  cavity  resonant 
frequency,  loss  due  to  diffraction  and  other  loss  mechanisms,  and 
spatial  distribution  of  the  transverse  mode  l  ,  respectively,  and 
each  mode  has  taken  on  a  factor  pi  ,  the  atomic  dipole  moment,  for 
simplicity.  D is  the  time  and  space  dependent  atomic  inversion  in 
tkv  cavity  and  D^.  is  the  steady  state,  space  independent  inversion  in 
the  absence  of  any  Held.  The  operators. 


2  dt  +  F1  +  r2 


P  \W>  =  r+  1<V  t 


operate  on  all  variables  to  their  right.  The  original  equations  included 
another  equation  for  the  polarisation  in  the  medium,  and  the  operators 
(3)  appear  in  (1)  and  (2)  because  the  polarization  has  been  eliminated 
from  the  equations.  In  (3),  and  are  the  decay  constants  for  the 
upper  and  lower  atojnic  or  molecular  states  a  and  b  in  Figure  1,  to 

3l 

is  the  resonant  frequency  and  T  is  the  decay  constant  for  the  polariza¬ 
tion.  Ir.  (1)  and  (2),  the  sums  m  and  n  are  over  ail  oscillating  inodes 
and  the  cum  M  is  over  all  the  atoms  or  molecules  in  the  cavity  and 
will  ultimately  be  replaced  by  a  volume  integral. 

We  make  the  substitutions  (4)  into  (i)  and  (2), 

b'  =  b.e'lV/ 

X  X 


r.  i&At 


where  v^  is  now  some  frequency  other  than  the  cold  cavity  frequency 

of  mode  f.,  and  b^.  its  slowly  varying  mode  amplitude.  Here,  A  -  \^-Vj 

and  'D^  depends  only  on  position.  In  the  equations  resulting  from 

substitution  of  (4)  into  (1)  and  (2),  wo  match  time  harmonics  and  obtain 

* 

the  equations  of  motion  for  the  amplitudes  of  the  two  oscillating  modes 
and  the  amplitudes  of  the  harmonics  of  the  inversion,  (all  other  modes 
have  zero  amplitude. ) 

*1  *■  t‘4,  -  f  1  *1  -  ppKTj  r  f  +  u2uiViM1  <5' 
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where, 


and 
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i  -■=  1,2. 


Tne  set  (7)  -  (9)  may  be  solved  by  perturbation  theory  assuming, 

D  >  D1^  >  D*yl>  Tf* 
w  o  ±1  ±2,  etc. 


and  v/n  obtain, 
.M 
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-dw  u2  vx  bz  bl  A21ai 
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Then  (13)  In  (5)  yields, 


I  Yj  j  /DjOjI  ^ 

J  ' v bi  i  1  ■  'r  ■  +iu^Fi^ i  2a  2(o-r / 


^  A21(n  _  ^ 

P{li)  \n^!u1l2{b1!571)(0)i  ju2?^22(o,  fi+  (u.  |2 lb, irA?I (i)+|u2 f •  jb2  ^a2Z<D3j 

,  (15) 

For  the  amplitude  of  mode  1  and  we  obtain  a  similar  expression  fer  the 

implitudp  of  mode  2  by  substitution  of  (13)  into  (6).  (  ) means  average  over 

-he  entire  laser  volume.  Hie  steady  state  number  of  photons  in  mode  l 

is  then  given  by 

Y>  =E,  li  /DwKi2  \  } 

2  LP<  '  "-!  \(l * |Uj i"  |i); !^A ,;( 3 H I U2  !2 lb2 A2 (0 1)  /  j 


!>>/  nJOghiiif  /..t<n 
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i2 1 ,  I?, , 


.v’ith  a  similar  e^prtsfion  for  the  steady  state  number  of  photons  in  mode 
i.  Similar  expressions  to  (15)  and  (16)  may  be  written  for  the  frequencies 
j{  two  modes,  except  that  they  involve  fhc  imaginary  part  of  the  same 
-repression. 

In  ord"r  to  com  pa  ve  the  steady  state  number  of  photons  in  the 
two  modes,  it  will  he  necessary  to  evaluate  4he  spatial  integrations 
■ver  Inc  modes  and  U2  as  denoted  by  the  brackets  (  ).  T hi >  is 


*  ' 


where  we  ran  into  insurr  mount  able  cliff.'  unities  in  evaluating  the 

complicated  expressions.  Some  remarks  about  tb°  "results”  (15)  and 

{16)  are  in  order  however. 

2  2  ' 

3)  If  jbjj  and  jb^j  are  taken  to  be  approximately  constant,  and  the 
spatial  integrations  in  (15)  taken  to  be  a  constant,  then  (15)  reduces  to, 

v2 

bj  +  (A-  --  +  i&j)  bj=0  (17) 

where  *h*»  definition  of  A  is  obvious  but  complicated.  Equation  (17)  has 


the  solution, 


bj=bo  exp  (Re  A  -  i  i(Aj  +  Im  A)t  | 


Now,  will  be  stable  (i.  e.  not  grow  or  decay  exponentially)  if, 

=  ReA 

This  little  argument  is  evidently  the  source  . of  (16)  uses  the 
reasoning  of  (LLM). 

2)  The  saturation  terms  come  in  as  expected.  In  (7_.LM),  for  a  single 
mode, 


i  +  b|b2f 


and  here  we  have, 


i  a  term  like  - 


i  +  »  vibj 


another  term 
like 


[HUjf  Ibjl" Au(0)+|ti2  f  |b2  f  A 2?_(0»- 

! 


i>  .  i..  i2  i,  i2  ,  , 


{!■»  jUj  \c  |bi  f  An(l)4- \VZ  l2 1  b2 1 2 A  22  (1)] 


Term  ia)  is  like  a  total  saturation  term  and  (b)  is  like  a  cross 

saturation  term  for  the  effect  of  mode  Z  on  mode  1. 

3)  Is  seems  strange  to  us  that  the  steady  state  number  of  photons  in 

T*  2 

mode  1  is  defined  in  terms  of  the  intensities  jbjj  and  [bgl  of  *»oth 


modes. 


Fig.  1.  Simple  I-aser  Atomic  or  Molecular  Model 
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ABSTRACT 

We  show  that  the  criterion  for  single  transverse  mode  laser 
operation  can  be  reduced  to  the  criterion  that  the  intensity  of  the  fun¬ 
damental  mode  be  greater  than  some  constant  in  the  entire  laser 
volume.  We  also  obtain  amplitude  and  phase  distributions  fOT  low  loss, 
edge  rimmed  unstable  resonators  for  Fresnel  number  up  to  eight. 

The  resonators  with  the  most  uniform  amplitude  distributions  are 
the  most  likely  to  suppress  higher  order  transverse  mode#. 


TRANSVERSE  MODE  SUPPRESSION  USING 
RIMMED  UNSTABLE  RESONATORS 


I.  Introduction 

There  has  been  much  interest  in  obtaining  a  high- poster,  high 

Fresnel  number  laser  that  will  operate  in  the  single  lowest  order  transverse 

mode.  In  general,  lasers  will  operate  in  numerous  transverse  modes 

simultaneously,  or  they  will  oscillate  in  a  high  order  transverse  mode. 

There  have  been  a  number  of  attempts  to  force  the  laser  to  oscillate  in  the 

1  2 

lowest  transverse  mode.  Birch  and  Skinner  and  Geusic  have  suggested 
the  cats -eye -pinhole  arrangement  that  spreads  out  the  volume  of  the  funda¬ 
mental  mode.  In  this  way,  the  fundamental  mode  saturates  out  more  of  the 

atoms  so  that  few  atoms  remain  to  excite  the  higher  order  modes.  At  high 

2 

powers,  however,  it  has  been  found  that  it  is  easy  to  destroy  the  pinhole. 

It  it  therefore  preferable  to  use  no  intra-cavity  optics  or  apertures. 

3 

Siegrnau  and  Arrathoon  have  studied  unstable  resonators  in  an 

attempt  to  spread  out  the  mode  volume  of  the  fundamental  mode  and  thu_- 

"discriminate"  against  higher  transverse  modes.  The  discrimination  is 

based  on  the  relative  losses  of  the  various  modes  and  Siegman  found  the 

diffraction  losses  of  even  the  fundamental  mode  to  be  quite  high.  Hence, 

it  is  possible  to  operate  in  a  single  transverse  mode  if  high  diffraction 

tosses  and  high  beam  divergence  can  be  tolerated. 

4 

Recently,  Lax,  Louisell,  and  Me  Knight  (LLM)  have  written  dovm 
the  coupled  equations  of  motion  for  the  mode  amplitudes  of  a  laser  oscillating 
in  number  of  •  '*"4  and  they  have  found  a  criterion  for  any 

13 


higher  order  mode  either  to  be  quenched  or  driven  into  oscillation  by  a  . 

i 

single  strongly  oscillating  lowest  order  mode.  They  have  evaluated  this 

criterion  for  a  confccal  resonator  and  fouad  that  nomeroos  higher  modes  t 

* 

would  be  driven  into  oscillation  by  the  single  lowest  order  mode. 

In  this  paper,  we  show  that  the  (LLM)  criterion  for  die  oscillation 

■ 

of  only  the  fundamental  mode  may  be  reduced  to  the  single  criterion  that 
the  intensity  distribution  of  the  fundamental  mode  be  greater  than  some 
constant  over  the  entire  laser  volume.  We  then  investigate  die  mode 
structure  of  unstable  resonators  with  rims  in  order  that  we  might  obtain  the 
low  diffraction  loss  mode  distribution  that  best  satisfies  this  criterion,  at 
lea£t  cn  the  mirror  surfaces.  The  unstable  part  of  the  resonator  spreads 
out  the  mode  volume,  while  the  runs  help  to  keep  the  energy  from  spilling 
out  of  the  resonator. 

It  is  found  that  by  varying  the  resonator  configuration  one  can  find 
favorable  intensity  distributions  at  fairly  low  loss  for  Fresnel  numbers  from  Z 
to  8.  6.  It  i3  also  found  that  unless  one  is  rather  careful,  he  can  easily 
force  interior  zeros  in  the  intensity  distribution  across  the  mirrors.  The 
appeararce  of  zeros  in  the  mirror  intensity  distribution  becomes  more 
sensitive  to  changes  in  mirror  configuration  at  higher  Fresnel  number. 

II.  Special  Case  of  (LLM)  and  the  Rimmed  Unstable  Resonators 

The  analysis  of  (LLM)  shows  that  a  single  oscillating  fundamental 
mode  will  quench  any  other  Iran  eve.****  r.iudc  if  tb'-  quotient  Q  satisfies 


where  y^,  y ^  and  ere  the  looses  and  spatial  distributions  of  the 

strongly  oscillating  fundamental  mode  E  and  any  other  higher  order  mode  1, 
respectively.  D  ia  the  ptnap  distribution  end 

*  K(0)  =  1  ♦S|uE|2 

i 

where  S  is  a  saturation  factor  that  depends  on  die  intensity  of  mode  E  and 
how  close  its  frequency  is  to  atomic  resonance 

2ii2r  lbEi2 

r(o)  [  r2  +  (Vg-o^2] 


Here,  \i  is  the  atomic  dipole  moment,  b  is  die  amplitude  of  mode  E  and 
v_  is  its  frequency,  o  is  the  atomic  resonance  frequency  and  r(o)  is  a 
combination  ol  atomic  decay  constants.  <  >  means  integration  over  the 
volume  of  the  active  medium.  (LLM)  have  evaluated  (1)  for  a  confocal 
resonator  and  have  found  that  mode  E  would  cause  numerous  higher  order 
transverse  mrdes  to  break  into  oscillation. 


Without  specifying  the  functional  dependences  of  u^.  and  on  r,  it 
-is  possible  to  reduce  (1)  to  a  simple  expression.  If  M  =  max  |u  j2  and 

1L 

C  -  min  Ju  |",  then  (1)  becomes 


“3‘15 


£ 


where  we  her:  taken  y  -  y  and  D  to  be  independent  of  petition.  Now, 

£  C  w 

(Z)  will  be  satisfied  if 

C  «  min  |  «  P  >  some  constant  (3) 

E 

over  the  entire  laser  volume. 

Equation  (3)  is  a  quantitative  statement  of  the  notion  that  one  most 
spread  out  the  mode  volume  of  the  single  oscillating  fundamental  mode  in 
order  that  it  saturate  die  whole  laser  medium  and  leave  no  gain  for  die  other 
modes. 

In  order  to  satisfy  die  criterion  Pi,  at  least  over  the  mirror s,  and  also 
keep  the  diffraction  losses  to  a  minimum,  *e  have  investigated  numerically 
the  eigenvalues  and  eigenmodes  of  the  type  of  resonator  shown  in  Figure  1. 

V be  amplitude  distribution  that  best  satisfies  (3)  and  keeps  the  diffraction 
losses  to  a  minima.:!  will  be  the  one  that  is  most  uniform  across  most  of  the 
mirror  and  then  falls  rapidly  to  zero  near  the  mirror  edge. 


JIT,  Integral  Equations 

The  resonator  in  Figure  1  is  equivalent  to  an  infinite  lens  system,  a 

part  r?  which  is  shown  in  Figure  2.  The  complex  amplitude  distributions 

n  (r)  are  defined  in  the  plane  immediately  preceeding  each  lens.  The 
n 

distribution  in  front  of  lens  n+1  that  arises  from  the  distribution  »u  'rent  of 
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lens  r,  in  the  case  where  d  »  jr|,  is, ' 


•where  the  integration  is  over  tibe  surface  preceding  lens  n.  We  make  use 
of  Use  follo»ing  to  reduce  (4)  to  specific  form: 

1}  The  phase  shift  suffered  by  any  ray  in  traversing  Use  distance  d,  in  a 
cylindrical!?  symetrre  system,  in  the  Fresnel  approximation,  is 


)  ~  Jd  +  (r’coa  r.*cos  4)  *"  (r*  sin^1  -  rsin^}2 


2.  2 


r  i  + r  *  r~  ~ 2r** cos  -  4*)  i 

L  ,.2  j* 


2)  The  phase  shift  associated  with  traversing  lens  n  is. 


kO,(r’)  =  — 
L  b‘ 


r  <  r 
—  o 


v>  *  k(d-4-; 


r'  >  r 
—  o 


where  a  must  be  adjusted  so  that 


o  L  b  ,  ,  J 


3)  To  form  an  integral  equation,  we  let  u  ,  =  vu  (r)  where  y  is  a 

n+I  u 

lied. 

complex  number  that  absorbs  all  uni-modual  phase  factors  like  e  and  i. 
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v  is  related  *o  the  diffraction  loss  of  the  leu*  system  per  lens  tra«*i  i. 


L  *  J-ivf- 


1)  Vf*z  stipulate  that  we  are  only  interested  in  modes  that  are  independent  of 
4-.  This  permits  os  to  do  the  4'  integration  with  the  usual  Bessel  function  as 
the  result. 


Y«*(r; 


j£q  cation  (4)  becomes  the  integral  equation, 

I  •*  ^■frZ+(]+2d/b')rl2] 

-  i  :  ufr’Je  '  J  J  f  t  rr^  r'  dr- 

°  (1 
k  ikd  frl  XdLr  +{1  "2d/b,r‘  J  -k 

+  -r  e  u(r')  e  J  {  —  rr*  )r‘  dr* 

a  J  o\d  / 


where  J  is  the  seroth  order  Bessel  function.  For  convenience  we  reduce 
o 

x 

(?)  t  ?.oro»?.Iiaed  form.  Letting  x  =  —  ,  we  obtain, 

T1 


i  \  ■>  f  f*  /  «  t  /)  mi  i  itrNl[x2+  (l+Zd/b‘)y2  3 

yu(x)  -  iuM  ^  j  ufy)  Jo(2trNlxy)e  7  y*jy 


ZittNof  J -  —  ^  1  inNlfx2+(l  ~2~)y2  i 

,  \b  b’  yr  .  .  L  V  b  J7  j 

+  e  I  u(y)  Jo(2trNlxy)e  yayj 


ro  n^1  To  ri 

where  J  -  —  =  /  and  NO  =  — —  and  N1  =  T  are  the  Fresnel 

r.  J  >d  Xd 

t  v 


numbers  ox  the  inner  and  cuter  mirrors,  respectively. 


In  terms  of  kernels,  (8)  is 


IS 


-6- 


SL 


i 


r 


q  * 

Yu(x)  =  J  K,(x,y)  u(y)dy  -r  J  ^(x.y)  o{y)  dy 
©  '  *~ 


(91 


where 


Kjfcc.y)  = 


H,Nl[x2+(l+|i)y2] 


2eKlJo(2i^ixjr)  e 


(10) 


J^{x;y}  -  2nNlJ^(2tiNixy>  e 


mm[*2+(l-f£)y2] 


T  * 


IV.  Method  o:  Solution 


¥/e  solve  the  integral  equation  {9)  by  approximating  the  integrals  on 

6 


the  right  band  side  of  the  equation  by  Legrendre -Gauss  quadrature.  There, 
we  represent  the  first  integral  by  a  sum  c £  M  terms  and  the  second  integral 
by  a  sum  of  tv  terms.  In  srdsj  to  use  this  type  of  quadrature,  we  must 
transform  the  ranges  of  integration  of  both  integrals  to  (-1, 1}.  The  details 
»r*'  rather  tedious  and  will  not  be  presented  hers.  The  integral  equation 
*9)  is  thereby  replaced  by  the  set  of  algebraic  equations. 


M 
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Yufaj) 


T>  K  (x.,x.)  B.  u(x.)  +  Tj  K  (x.,x.)  B.  u(x.) 
j-1  A  1  J  J  J  j=l  *  1  J  J  1 


(ID 


’vherc  tifrr j ;  rhe  complex  amplitude  at  unequally  spaced  mirror  positions 


k..  The  x.  arc  related  to  the  zeros  of  the  Legendre  polynomials  P  (xl  and 
j  J  M 

n 


P^(xj,  and  B;  are  :abu!ated  weight  functions.  *  Equation  (11)  is  equivalent 
to  the  matrix  equation,  19 
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9 


Kn  =  y  a 

where  K  is  ac  (M  x  (M+N)  complex  matrix  «hese  elements  are 

K  B.  for  j<  M 

*  1  J  J 

for  j  >  M 

and  the  vector  a  is  the  complex  amplitude  evaluated  at  the  various  x  .  The 

m!rir  elas-mb  are  easily  generated  on  computer  and  we  then  use  the  IBM 

SHAPE  subroutine  ALLMAT  to  find  the  eigenvalues  and  right  eigenvectors 
~  8 

of  the  matrix  K.  The  right  eigenvectors  are  then  ice  eigenmodes  cf  the 
cavity. 

Asa  check  on  the  mathematical  formulation  and  the  computer 
program,  we  have  plotted  the  results  of  the  special  cases  of  plane  parallel 
(d/b  =  d/b’  -  0}  and  esnfocal  (d/b  -  1,  d/b*  =  -1)  resonators  and  compared 
then*,  with  the  results  of  the  wave  launching  method  of  Fox  and  Li  in 
Figure  3.  Evidently,  the  agreement  with  the  old  results  is  quite  good 
and  *>e  *.n*;  believe  the  results  for  all  reasonable  values  of  d/b  and  d/b'. 

V,  Results  and  Discussion 

Having  convinced  ourselves  that  the  computer  program  works  for 
widely  varied  d'b  and  d/b1,  we  can  now  proceed  to  investigate  the  behavior 
cf  the  irirrur  amplitude  distributions  as  a  function  of  d/b  and  d/b’.  In 

in  fig,;rr  3  ac  ha  shown  only  the  amplitude  distributions,  but  the  phase 
attributions  compare  equally  well. 

&) 


order  to  get  some  qualitative  feel  for  hoe  «e  expect  die  distributions  to 
depend  on  she  parameters,  it  is  necessary  to  solve  (8)  in  some  approximate 
sense.  Unfortunately,  cot  ody  is  die  integral  ec nation  (8)  composed  cC  two 
parts,  »he  parts  are  ucn>hermitian,  oon-symetric  and  non-degenerate. 
Furthermore,  sir.ee  we  are  interested  in  large  as  well  as  small  Fresnel 
numbers,  it  is  pointless  to  expand  J^ttNlry)  in  a  power  series  in  its 
argument  (the  usual  trick  for  obtaining  a  degenerate  kernel;!^  We  have 
spent  considerable  t  me  trying  to  obtain  any  kind  of  approximate  analytical 
solution  to  (3)  and  hive  been  unsuccessful.  Therefore,  any  qualitative  feel 
for  ho%  the  distriira'doos  will  depend  cc  the  parameters  is  completely 
missing  and  we  must  proceed  on  a  trial  and  error  basis.  Hence,  it  will  be 
rattier  difficult  to  explain  a  logical  progression  to  the  desired  results,  and 
we  will  only  be  able  to  indicate  how  we  obtained  the  results  and  hope  to 
show  rfceir  apparent  importance.  We  begin  with  low  Fresnel  number. 

A.  Low  Fresnel  Number 

lit  Figure  4  we  saw  a  series  of  mirror  amplitude  distributions  for. 

1  =  2  and  MO  =  .25.  (The  intersection  of  the  interior  and  exterior  mirrors 
is  indicated  by  the  hash  marks  on  the  abscissas  on  all  of  the  remaining 
fig-_re5.  The  diffraction  Iocs  of  the  lowest  loss  mode  is  L  and  that  of  the 
next  lowest  loss  n  ode  is  occasionally  listed  in  parenthesis.  }  This  resonator 
is  definitely  of  cor.f oral- type  with  an  interior  spherical  hump.  We  held  the 
"ccnfocai  paiax^eti- r:i  d/b  constant  at  .2  and  steadily  Increase  d/b*  to  force 
rr-ore  of  the  energy  to  the  outside  of  the  reconatoi.  The  procedure  seems 
to  be  v/otking  relatively  well  until  an  interior  zero  in  the  amplitude  distribution 


*_ 

« 

begins  to  form  and  at  d/b*  -  4.  die  distribution  becomes  quite  unacceptable. 

To  remedy  tins  situation,  *e  went  to  fine  ana  table-type  resonator 
with  rims  shown  in  Figure  5  for  N1  =  2.,  NO  =  1.  5.  Here,  die  first  attempt 
yielded  an  acceptable  solution  but  the  diffraction  loss  was  rather  high  {25%/ 
pass).  Therefore,  we  increased  d/b  (made  die  rims  more  carved)  to  try 
to  decrease  the  diffraction  loss.  At  d/b  =  .4,  d/b*  -  •  1  we  find  a  rather 
acceptable  solution  with  a  loss  of  17%  per  pass.  It  ie  interesting  to  note 
that  the  amplitude  ripples  are  of  higher  frequency  and  of  lower  amplitude 
for  NO  HI  (unstable  with  rims)  and  cf  lower  frequency  and  higher  ampli¬ 
tude  for  NO  «  N1  (coafocal  with  interior  bumps).  The  unstable  with  rims 
are  therefore  of  more  interest  to  us  because  large  amplitude  variations 
across  the  mirror  will  undoubtedly  cause  zeros  and  the  distribution  will  not 
satisfy  th3  basic  criterion  (3).  On  the  other  baud,  rapid  amplitude  variations 
ars  acceptable,  provided  they  are  of  low  rned elation. 

The  phase  distributions  in  either  case  did  not  change  appreciably 
fvcm  nai  to  run  so  *.v«r  have  plotted  only  a  typical  phase  distribution  in  the 
lower  right  hand  comers  of  Figures  4  and  5.  Neither  phase  distribution 
varies  appreciably  over  the  entire  mirror  surface,  but  we  do  note  that  the 
phase  /arios  mere  rapidly  neir  the  origin  for  Nt  -  Wl. 

B.  Medium-  Fresnel  Number 

We  found  similar  behavior  for  Nl  =  4.  as  for  Nl  =  2.  Namely,  it 
was  easier  to  find  accspabie  distributions  for  NO  -  Nl,  shown  in  the  left, 
column  of  Figure  7,  rather  than  for  NO  <<  Nl,  shown  in  Figure  6.  We 
also  found  that  for  NC  and  Nl  integer,  the  amplitude  distribution  could 


'"-  — -.t-L—x-.t  -. 


have  rather  abrupt  change*  for  email  changes  in  d/b  or  d/b’.  We  feel  that 
this  high  sensitivity  to  d/b  end  d/b1  is  probably  due  to  die  fact  that  for  small 


d/b'  die  exponentials 


and 


i«Nl(l  +  §£)y2 

*«0  -  r  >2 


in  ibc  kernels  in  equation  (8)  oscillate  oyer  very  nearly  an  integral  number 
of  half  cycles  over  the  ranges  (0,q)  and  (q,  1).  We  therefore  felt  that  we 
would  have  more  success  for  NO  and  Nl  not  integers  and  this  is  borne  out 
by  die  distributions  shown  in  the  right  hand  column  of  Figure  7.  We  feel 
that  the  second  distribution  in  this  column  will  be  successful  in  suppressing 
higher  order  transverse  modes  at  this  Fresnel  number.  Note  that  the 
diffraction  loss  is  only  4.  5%  per  pass. 

C.  Higher  Fresnel  Number 

In  Figure  8  we  have  plotted  the  amplitude  distributions  for  rather 
high  Fresnel  number,  and  we  feel  that  the  distribution  at  the  bottom  of  the 
left  column  will  be  successful  in  suppressing  higher  order  modes.  Hera, 
the  distributions  become  very  sensitive  to  changes  in  d/b  and  d/b'  and  for 
typical  lasers  a  favorable  distribution  would  probably  require  that  the 
miners  have  radii  of  curvatures  in  the  hundreds  of  meters.  We  feel, 
however;  that  provided  the  mirror  is  spherical  to  1/20  *h  of  a  wave,  that 
the  distribution  will  be  approximately  the  same  even  for  large  tolerances 
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on  the  radius  (L  e. ,  400  or  450  meters). 
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VI.  Conclusions 


Although  the  unstable  resonators  of  Siegman  3nd  Arrathoon  can  be 
shown  to  have  high  discrimination  against  higher  order  transverse  modes, 
this  discrimination  does  not  guarantee  suppression  of  these  higher  modes. 
Indeed,  as  shown  in  {LLM),  if  the  gain  for  a  particular  mode,  caused  by 
the  presence  of  a  strongly  oscillating  fundamental  mode,  exceeds  the  loss, 
the  higher  mode  will  break  into  oscillation,  independent  of  the  relative  los3e3 
of  the  two  modes.  The  analysis  of  (LLM),  has  provided  us  with  a  criterion 
for  the  suppression  of  any  higher  order  mode  by  the  fundamental  mode,  all 
gains  included. 

In  this  paper,  we  have  (1)  shown  that  this  criterion  can  be  reduced  to 
a  quantitative  statement  of  the  notion  that  one  needs  to  spread  out  the  mode 
volume  of  the  fundamental  mode  in  order  to  suppress  higher  order  modes, 
and  (2)  solved  numerically  for  the  low  diffraction  loss  eigenmodes  of 
rimmed  unstable  resonators  in  order  that  we  might  find  the  mode 
amplitude  distributions  that  best  satisfy  this  criterion  over  the  mirror 


surfaces. 


24 


References 


3.  J<  M.  Burch,  “Design  of  Resonators",  in  Quantum  Electronics,  ed. 

P.  Grivet  and  N.  Bloembergen  (Columbia  University  Press,  New  York, 
1964),  Vol.  2,  pp.  1187-1202. 

2.  J.  G.  Skinner  and  3.  E.  Geusic,  "A  Diffraction  Limited  Oscillator",  in 
Quantum  Electronics,  ed.  P.  Giivetand  N.  Bloembergen  (Columbia 
University  Press,  New  York,  3  964),  Vol.  2,  pp.  1437-1444. 

3.  A.  E.  Siegman  and  R.  Arrathoon,  "Modes  in  Unstable  Optical 
Resonators  and  Lens  Waveguides",  IEEE  J.  Quantum  Electronics, 

Vol.  3,  pp.  156-163,  April  1967. 

4.  M.  Lax.  W.  H.  Loiusell  and  V7.  B.  McKnight,  to  be  published  in  IEEE 
J.  Qnanti.m  Electronics. 

5.  M,  Bor?,  and  E.  Wolf,  Principles  of  Optics,  (Pergamon  Press, 

Now  York,  J  970),  Ch.  8. 

6.  F.  B.  Hilicbrand,  Introduction  to  Numerical  Analysis  (McGraw-Hill, 
New  York,  i%5),  Ch.  3. 

7.  C.  H.  Love,  "Abscissas  ans  Weights  for  Gaussian  Quadrature  for 
N  -•  2  (u  100  and  N  =  12?,  1  50,  175  and  200",  National  Bureau  of 
Standards  Monograph  98  (1966). 

8.  1L  L.  Sanderson  and  W.  Streifer,  "Comparison  of  Laser  Mode 
Calculations",  Applied  Optics,  Vol.  8,  pp.  131-136,  Jan.  196V. 

9.  A.  G.  Fox  and  T.  Li,  "Resonant  Modes  in  a  Maser  Interferometer", 

Bell  Sys.  Tech.  J. ,  Vol.  40,  pp.  453-488,  March  1961. 

10.  D.  E.  McCumber,  "Eigenmodes  of  a  Symmetric  Cylindrical  Confocal 

Laser  Resonator  and  Their  Perturbation  by  Output-Coupling  Apertures", 
Bell  Sys.  Tech.  J. ,  Vol.  44,  pp.  333-363,  Feb.  1965. 

25 


-13- 


FIGURE  CAPTIONS 


Figure  i 
Figure  2 

Figure  3 

Figure  4 

Figure  5 
Figure  4 
Figure  7 

Figure  8 


Unstable  resonator  with  rims 

Lens  system  that  is  equivalent  to  j  esonator  shown  in 
Figure  1 


Comparison  with  Fox  and  Li  wave  launching  method  for  plane 
parallel  and  confocal  resonators. 

— - - - - —  —  Fox  and  Li 

xxxxxxxxx  Ours 

Mirror  amplitude  and  phase  distributions  for  varicus  d/b, 
d/b'  at  low  Fresnel  number.  Confocal-type  with  Interior 
spherical  hump. 

Mirror  amplitude  and  phase  distributions  for  various  d/b 
and  d/fc1  at  low  Fresnel  number.  Unstable  with  rims. 

Mirror  amplitude  and  phase  distributions  at  medium  Fresnel 
number.  Confocal  with  interior  spherical  humps. 

Amplitude  distributions  at  medium  Fresnel  number  for  integer 
and  non-integer  Fresnel  numbers.  Unstable  resonator  with 
rims. 


High  Fresnel  number  amplitude  and  phase  distributions. 
Unstable  with  rims. 
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